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Abstract 

Doping and disorder are inseparable in the superconducting cuprates. Assuming the 
simplest possible disordered doping, we construct a semiphenomenological model 
and analyze its experimental consequences. Among the affected experimental quan- 
tities are the ARPES spectra and thermodynamic properties. From our model we 
make a prediction for the width of the local superconducting gap distribution with 
the only experimentally unknown parameter being the superconducting correlation 
length. Thus, our model provides a direct way of determining the superconducting 
correlation length from a known experimental gap distribution. 
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1 Introduction 

Ground state of all cuprate superconductors is a Mott insulator at half filling. 
To render this state a superconductor one has to dope it with carriers (holes for 
p-type superconductors such as LSCO and Bi2212, and electrons for n-type, 
e.g NdSCO). Below we will focus on hole doped superconductors, although 
discussion about intrinsic disorder is relevant for electron doped materials as 
well. In the low doping limit, dopants, such as Sr in LSCO and oxygen in 
Bi2212 are certainly randomly distributed across the sample without form- 
ing any regular array. The central point we argue for in this paper is that 
effect of carriers and intrinsic disorder produced by the very same doping are 
inseparable in underdoped and possibly overdoped cuprates. 
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The microscopic inhomogeneity of the superconducting state state has impor- 
tant consequences for the average properties obtained in the bulk measure- 
ments. Immediate conclusion we draw is that even with the perfectly grown 
crystals doping will produce artificial impurities that disorder the system. The 
affected properties include the average tunneling spectra, as the ones obtained 
by ARPES and measured specific heat. Another implication is that since the 
state is inhomogeneous, there should be no well defined superconducting tran- 
sition temperature, much like in the conventional granular superconductors [1]. 
This "granularity" should also clearly enhance the electric and the thermal 
transport above T c . 



2 Density of States and gap distribution in the presence of doping 
induced disorder 



To analyze the effects of doping-induced disorder, we use a minimal model 
of randomly distributed dopants in the reservoir layer. The superconducting 
order parameter can only change on the length scale of the superconducting 
coherence length, £. Hence we can coarse-grain the Cu-0 plane into patches 
of the size £ and on these patches define the effective dopant concentration as 

, s Ni(r) 

nc(r) = -p-A (1) 

where iVj(r) is the number of dopants in the adjacent reservoir plane which 
happen to be within a circle of radius £ around the point r (Figure 1). 




Fig. 1. The effective local doping at a particular point is defined through the random 
number of dopants which happen to be in the ^-vicinity of this point. 

Assuming that the positions of the dopants are uncorrelated, which is expected 
to hold particularly at low doping levels, the numbers N% obey the Poisson 
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distribution, with a probability density, 



P(iV) = ^|^exp(-< 2 ). (2) 



Here, the average dopant density, n, is a function of the nominal doping frac- 
tion, x = abn, with a and b the in-plane lattice constant. An example of 
Poisson distribution is shown in Figure 2. 
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Fig. 2. Poisson distribution with the average equal to 3. The distribution is defined 
only on non- negative integers. 

The main characteristic of the Poisson distribution is that its dispersion is 
equal to its average value. In the limit of large average, N = nC, 2 ^> 1, the 
Poisson distribution reduces to a Gaussian. In this limit, the probability dis- 
tribution of the effective local doping fractions becomes 
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The distribution of the local doping should lead to a corresponding distribution 
of the gap values over the regions of size £. Assuming that the distributions 
of doping is sufficiently narrow, the value of the local gap as a function of the 
local doping can be linearized as 

A(x) = A - Aix. (4) 



To the first order approximation, the parameters Ao and Ai can be extracted, 
e.g. from the ARPES data which determines the average gap values as a func- 
tion of the average doping. Then, from Equations (3) and (4), the probability 
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distribution for the local superconducting gaps follows: 
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(5) 



The probability distribution of the local values of the gap A(x) is completely 
specified up to the unknown parameter £. Hence, one of the applications of the 
local gap statistic, and in particular the distribution width is the determination 
of the doping-dependent correlation length, £. 

To illustrate the smoothing effect of the gap fluctuations on the bulk properties 
we compute the average density of states (DOS), as measured by photoemis- 
sion. In a conventional <i-wave BCS superconductor with the amplitude of the 
superconducting gap A, the density of state is given by 



Here, N is the electronic DOS at the Fermi level in the normal state and 
K(x) is the complete elliptic integral of the first kind. 

To determine the average DOS is the presence of the gap amplitude fluctua- 
tions, we compute the integral 



The average DOS, N(u), clearly preserves the some rule / dujN{uj) = 1 due 
to normalization of -P(A). 

In Figure 3 we show the DOS curves for various widths of the gap distribution. 
In the "clean" case one indeed observes a divergence in DOS at the energy 
equal to the amplitude of the gap. For a finite-width gap distributions, the 
divergence is replaced by a hump near the average gap value. For wide gap 
distributions, such that the distribution width is comparable to the magnitude 
of the gap itself (case a — 1 in Fig. 3), normal regions can be expected to 
appear in the sample. Such regions give a finite contribution to the DOS at 
the Fermi level. This contribution is similar in spirit to the well-known effect 
of non-magnetic impurities on the DOS in superconductors [2,3]. 

For the optimally-doped Bi2212, the observed width of the gap distribution 
is about 10% of the average [4]. As is clear from the Figure 3, such width 
corresponds to a rather smooth peak, much like the one observed in ARPES 
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[5]. 
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Fig. 3. The average BCS <i-wave density of states (DOS) for different widths of 
the gap distribution, a. Such spatial averaging may explain the smooth features 
commonly observed in the photoemission of optimally and underdoped cuprates. 
The finite DOS at the Fermi level for large a is caused by the zero-gap (normal) 
regions. 
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